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In general, the field equations of /(i?) theory coupled to a matter field are very complicated 
and hence it is not easy to find exact analytical solutions. However, if one considers traceless 
energy-momentum tensor for the matter source as well as constant scalar curvature, one can derive 
some exact analytical solutions from f{R) theory coupled to a matter field. In this paper, by 
assuming constant curvature scalar, we construct a class of charged rotating black string solutions 
in /(7?)-Maxwell theory. We study the physical properties and obtain the conserved quantities of 
the solutions. The conserved and thermodynamic quantities computed here depend on function 
f'{Ro) and differ completely from those of Einstein theory in AdS spaces. Besides, unlike Einstein 
gravity, the entropy does not obey the area law. We also investigate the validity of the first law 
of thermodynamics as well as the stability analysis in the canonical ensemble, and show that the 
black string solutions are always thermodynamically stable in /(i?)-Maxwell theory with constant 
curvature scalar. Finally, we extend the study to the case where the Ricci scalar is not a constant 
and in particular R — R{r). In this case, by using the Lagrangian multipliers method, we derive an 
analytical black string solution from f{R) gravity and reconstructed the function R{r). We find that 
this class of solutions has an additional logaritmic term in the metric function which incorporates 
the effect of the f{R) theory in the solutions. 
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I. INTRODUCTION 

There has been considerable attentions in the past years in modified gravity theories, specially f{R) theory which 
is one of the encouraging candidates for explaining the current accelerating of the universe expansion [l|, y] (see also 
Q for a comprehensive review on f{R) theories). In fact /(i?) theories can be regarded as the simplest extension of 
general relativity. Many f(R) models have passed all the available experimental tests and fit the cosmological data. 
To prevent a ghost state, f'{R) > for R > Rg is required d, Q- f"{R) > for i? > i?o, is needed to avoid the 
negative mass squared of a scalar-field degree of freedom (tachyon) [3|. f{R) -^ R — 2A iov R > Rq, is required for 
the presence of the matter era and for consistency with local gravity constraints Q. It was shown that f{R) theories 
can be considered as general relativity with an additional scalar field that provide new insight in the two cases of 
Brans-Dicke theory with wq = and ujq = —3/2 [y]. 

There have been a lot of works in the literature attempting to construct static and stationary black hole solutions 
in f{R) gravity theories. One may expect that some signatures of black holes in f{R) theories will be in disagreement 
with the expected physical results of Einstein's gravity. In [3| the authors studied general solutions in f{R) theory 
using a perturbation approach around the Einstein-Hilbert action. In § black hole solutions were found by adding 
dynamical vector and tensor degrees of freedom to the Einstein-Hilbert action. Also, the transition from neutron 
stars to a strong scalar- field state in f{R) gravity has been studied in Q. Physical properties of the matter forming 
an accretion disk in the spherically symmetric background in f{R) theories were explored in [l3|. In Ref. [H] the 
construction of traversable wormhole geometries was discussed in f{R) gravity. The Schwarzschild-de Sitter black 
hole like solutions of f{R) gravity were obtained for a positively constant and a non-constant curvature scalar in 
[13 | and [lj|, respectively. A black hole solution was obtained from f{R) theories by requiring the negative constant 
curvature scalar Q. If 1 + /'(Ro) > 0, this black hole is similar to the Schwarzschild-AdS (SAdS) black hole. It was 
argued that f{R) and SAdS black holes have no big difference in thermodynamic quantities when using the Euclidean 
action approach and replacing the Newtonian constant G by GcS = G/{1 + /'(i?o)) \f\- I^ i^ ^-l^o interesting to study 
black hole solutions in f{R) theory coupled to a matter field. In general, the field equations of f{R) theory coupled to 
the matter field are very complicated and hence it is not easy to find exact analytical solutions. In order to construct 
the constant curvature scalar black hole solutions from f{R) gravity coupled to the matter, the trace of its energy- 
momentum tensor T^^ should be zero [lj|. Two examples for the traceless T^j^ are Maxwell and Yang-Mills fields 
which were studied in J14l |l5| . Thermodynamics and properties of these solutions were also studied in ample details 
[14] . It was found that these solutions are similar to the Reissner- Nordstrom- AdS (RNAdS) black hole when making 
appropriate replacements jl4| . The Kerr- Newman black hole solutions with non-zero constant scalar curvature in 
/(i?)-Maxwell theory, their thermodynamics, as well as their local and global stability were also studied in [1(^]. 

In this paper we would like to continue the investigation on the f{R) black holes, by constructing a new class of 
charged rotating black string solutions in i? + /(_R)-Maxwell theory with constant curvature scalar. The traceless 
property of the energy-momentum tensor of the Maxwell field plays a crucial role in our derivation. With assumptions 
i?o < and 1 -I- /'(i?o) > our solution is similar to charged black string solution in AdS space with suitable replacing 
the parameters. We will also suggest the suitable counterterm which removes the divergences of the action. We calcu- 
late the conserved and thermodynamic quantities of these black strings by using the counterterm method. We obtain 
a Smarr-type formula for the mass of the black string and check the validity of the first law of thermodynamics. We 
perform the stability analysis in the canonical ensemble and show that the black strings are always thermodynamically 
stable in /(i?)-Maxwell theory with constant curvature scalar. Finally, we extend the study to the case where the 
Ricci scalar is not constant and in particular R — R{r) and derive an analytical black string solution. 

II. FIELD EQUATIONS AND SOLUTIONS 

We start from the four-dimensional R -I- f{R) theory coupled to the Maxwell field 

-^{R + f{R)-F^,Fn~^ [ d^xV^Q{h), (1) 

where R is the Ricci scalar curvature, F^i, = d^A^ — d^A^ is the electromagnetic field tensor, and A^ is the electro- 
magnetic potential. The last term in Eq. (jlj is the Gibbons-Hawking boundary term. It is required for the variational 
principle to be well-defined. The factor 9 represents the trace of the extrinsic curvature for the boundary dA4 and 
h is the induced metric on the boundary. The equations of motion can be obtained by varying the action ([I} with 
respect to the gravitational field 5^^ and the gauge field A^ which yields the following field equations 

R^, (1 + fiR)) - ig^,(i? + /(i?)) + (5^,V2 - V^V,) /'(R) = SttT^,, (2) 



V^F^" = 0, (3) 

with the energy-momentuni tensor 

T,. = ^ {f,,F,^ - ^g^^F,,F'^^ . (4) 

The above energy-momentuni tensor is traceless in four dimension, i. e., T^ = 0. As we mentioned aheady this 
property plays an important role in our derivation. In Eq. 1^ the "prime" denotes differentiation with respect to 
curvature scalar R. Assuming the constant curvature scalar R ~ Rq, the trace of Eq. © yields 

i?o(l + /'(i?o))-2(i?o + /(i?o)) = 0, (5) 

Solving the above equation for negative Rq, gives 

Substituting the above relation into Eq. ©, we obtain the following equation for Ricci tensor 

1 f fjRo) \ , 2 ^ 

'" " 2^^" [f'{Ro) -lj^l + f'{Ro)'"" ^ ' 

Now, we want to construct charged rotating black string solutions of the field equations ([2|) and ([3|) and investigate 
their properties. We are looking for the four-dimensional rotating solution with cylindrical or toroidal horizons. The 
metric which describes such a spacetime can be written in the following form [13, [l3 



ds^ = -Nir) {Edt - ad(f>f + r^ (^^dt - Edcj^ " + -^j-^ + ^dz"^ 



dr^ 

mr) ' P 



2' = l+p, (8) 

where a is the rotation parameter. The function N(r) should be determined and I has the dimension of length which 
is related to the constant Af by the relation P = — 3/Af. The two dimensional space, inconstant and r ^constant, 
can be (i) the flat torus model T^ with topology S^ x S^, and < (j> < 2tt, < z < 2ttI, (ii) the standard cylindrical 
model with topology R x S^, and < (f> < 27r, — cxd < z < cxd, and (iii) the infinite plane R^ with — oo < (j) < oo and 
— oo < z < oo. We will focus upon (i) and (ii). The Maxwell equation ([3]) can be integrated immediately to give 

F - '^^ 

i'tr — — , 
F^r = —^Ftr, (9) 

where q is the charge parameter of the black string. Substituting the Maxwell fields ([9|) as well as the metric ([U in 
the field equation ([2]) with constant curvature, the non-vanishing independent components of the field equations for 
a = reduce to 

(1 + f'iRo)) (2^^^ + 4.3^ + R.r^'^ - V = 0, (10) 

(1 + f'{Ro)) (4^"^^ + 4r2iV(r) + RqA + V = 0. (11) 

One can easily show that the above equations have the following solution 

N(r) = -^ + ^- ^r^ (12) 

^ ' r + (1 + /'(i?o))r2 12 "^ ' ^''> 

where m is an integration constant which is related to the mass of the string. One can also check that these 
solutions satisfy equations ©-(O in the rotating case where a 7^ 0. It is apparent that this spacetime is similar with 
asymptotically AdS black string. Indeed, with the following replacement 

^' ■ Q2 (13) 



(l + /'(^o)) 
Ro 



4 '^ (14) 
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FIG. 1: The function N{r) versus r for m — 2, f'{Ro) = 2 and q = 1. Ro = 12 (bold line) and Rq — —12 (continuous line). 
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FIG. 2: The function N{r) versus r for m = 2, g = 1 and Ro = -12. /'(-Ro) = 0.5 (bold hue), /'(i?o) = (continuous line) 
and f'{Ro) — —0.5 (dashed line). 

the solution reduces to the asymptotically AdS charged black string for A = —3/P [17i|. Next we study the physical 
properties of the solutions. The Kretschmann scalar for this solution is given by 



^^u\k^ 



u\k 



3r8(l + /'(i?o))2 



^"(tTT^oV^ + 18m2)(l + f{Ra)f - 36mrg2(i + f{R,,)) + 2lq^ 
Id 



(15) 



When r — > 0, the dominant term in the Kretschmann scalar is 56g^/[(l -I- /'(i?o))^''*]- Therefore we have an essential 
singularity located at r = 0. The Kretschmann scalar also approaches Rq^ /6 as r — > oo. As one can see from Eq. 
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FIG. 3: The function N{r) versus r for m = 2, f'{Ro) — —2 and q ~ I- Ro = 12 (bold line) and Ro ~ —12 (continuous line). 



(rT2)) . the solution is ill-defined for /'(i?o) = ^1- The cases with /'(i?o) > — 1 and f'{Ro) < — 1 should be considered 
separately. In the first case where /'(i?o) > ^1: there exist a cosmological horizon for i?o > 0, while there is no 
cosmological horizons ii Rq < (see fig. [T]). Indeed, for 1 + f'{Ro) > and i?o < the black string can have 
two inner and outer horizons provided the parameters of the solutions are chosen suitably (see fig. [2]) . In the latter 
case {/'(Rq) < —1), the signature of the spacetime changes and the conserved quantities such as mass and angular 
momenta become negative, as we will see in the next section, thus this is not a physical case and we rule it out from 
our consideration (see fig. |3]). 

III. CONSERVED AND THERMODYNAMIC QUANTITIES 

Next, we calculate the conserved quantities of the solutions by using the counterterm method inspired by (A)dS/CFT 
correspondence [19|. The spacetimes under consideration in this paper has zero curvature boundary, Rabcd{h) = 0, and 
therefore the counterterm for the stress energy tensor should be proportional to h"'^. We find the suitable counterterm 
which removes the divergences of the action in the form 




(16) 

IdM 

where Rq < 0. Having the total finite action I = Iq + let at hand, one can use the quasilocal definition to construct 
a divergence free stress-energy tensor [20| ■ Thus the finite stress-energy tensor can be written as 



rpab 

^ Stt 



Qab _ Q^ab _ J _^,^ab 



(17) 



The first two terms in Eq. (|17l) are the variation of the action ([1]) with respect to hab, and the last term is the variation 
of the boundary counterterm (J16p with respect to hab- To compute the conserved charges of the spacetime, one should 
choose a spacelike surface B in dAi with metric aij, and write the boundary metric in ADM (Arnowitt-Deser-Misner) 
form: 

habdx^dx" = -N^dt^ + Oij {dip' + V^dt) {dip^ + V^dt) , 

where the coordinates Lp^ are the angular variables parameterizing the hypersurface of constant r around the origin, 
and N and V are the lapse and shift functions respectively. When there is a Killing vector field ^ on the boundary, 
then the quasilocal conserved quantities associated with the stress tensors of Eq. p?)l can be written as 

Q(a = / d'x^Tabn'^e, (18) 

JB 

where a is the determinant of the metric Uij , ^ and n° are, respectively, the Killing vector field and the unit normal 
vector on the boundary B. The first Killing vector of the spacetime is f = d/dt^ and therefore its associated conserved 
charge of the string is the mass per unit volume. A simple calculation gives 

M = I d^x^Ta.n'^e = ^^^'g"/^"" [1 + /'(i?o)] ■ (19) 

The second conserved quantity is the angular momentum per unit volume associated with the rotational Killing 
vectors <; = d/d(f> which can be calculated as 



J = / d^xV^TabU^" = ^ [1 + /'(i?o)] . (20) 

Jb o"^ 

For a — (E. — 1), the angular momentum per unit volume vanishes, and therefore a is the rotational parameters of 
the spacetime. Next we calculate the entropy of the black string. Let us first give a brief discussion regarding the 
entropy of the black hole solutions in f{R) gravity. To this aim, we follow the arguments presented in |21[. If one 
use the Noether charge method for evaluating the entropy associated with black hole solutions in f{R) theory with 
constant curvature, one finds |12l | 

5=^/'(i?o), (21) 



where A = Anr"^ is the horizon area. As a result, in f{R) gravity, the entropy does not obey the area law and one 
obtains a modification of the " area law". Motivated by the above argument, for the rotating black string solution in 
R + f{R) gravity, we find the entropy per unit length of the string as 

S=^-^[l + nRo)]. (22) 

Then we obtain the temperature and angular velocity of the horizon by analytic continuation of the metric. Although 
our solution is not static, the Killing vector 

X = dt + nd^ (23) 

is the null generator of the event horizon where D, is the angular velocity of the outer horizon. The analytical 
continuation of the Lorentzian metric by t — > ir and a -^ ia yields the Euclidean section, whose regularity at r = r+ 
requires that we should identify t ^ T + /3+ and cj) ^ (/)+i/3+r2+ where /3+ and il+ are the inverse Hawking temperature 
and the angular velocity of the horizon. We find 

1 UN{r) \ [RorX{l + f'{Ra))+Aq^ ^ 



^ 47te[ dr )^^^^ m7TE[l + f'{Ro)]rl ' ^^^^ 

^ = ^, (25) 

where we have used equation N{r^) = for omitting the mass parameter m from temperature expression. Since 
1 + f'{Ro) > 0, therefore the temperature is non negative provided 

R,rl{l + f{Ro)) < -V -^Ro< -^AJ^^JT^]^^ (26) 

where the equality holds for extremal black string with zero temperature. The next quantity we are going to cal- 
culate is the electric charge of the string. To determine the electric field we should consider the projections of the 
electromagnetic field tensor on special hypersurface. The normal vectors to such hypersurface are 

TV' ' AT' ^ ' 

where A^ and V^ are the lapse function and shift vector. Then the electric field is E^ = g^PPpjyu'^, and the electric 
charge per unit length of the string can be found by calculating the fiux of the electric field at infinity, 

Q = , ^'^ (28) 

inl^l + f'{Ro) 

The electric potential U, measured at infinity with respect to the horizon, is defined by [22| 

U = A,x''\r^o.~ApX''l=r^, (29) 

where x is the null generator of the event horizon given in Eq. (j23p . One can easily obtain the electric potential as 

q 



U=^^l + f'iRo). (30) 

Then, we consider the first law of thermodynamics for the black string. In order to do this, we obtain the mass M as 
a function of extensive quantities S, J and Q. Using the expression for the mass, the angular momenta, the entropy 
and the charge given in Eqs. ([TO]) . (PU]). ([^ and (^5)) and the fact that N{r+) — 0, one can obtain a Smarr-type 
formula as 

MjS^J^Q)^ J( 3Z~1) ^3^^ 

where Z = S^ is the positive real root of the following equation: 



3VZ~1tt^Q^P[1 + /'(i?o)]^ - 2jTTJy/Z(l + f'{Ro))Sl + 3S^VZ~1 

= 0. (32) 



2Tr^VZ{l + f'{Ro))Sl 



One may then regard the parameters S, J and Q as a complete set of extensive parameters for the mass M{S, J, Q) 
and define the intensive parameters conjugate to S, J and Q. These quantities are the temperature, the angular 
velocities and the electric potential 






j[[7TQl{l + f'{Ro))f-3S^} 
3SI^Z{Z-1)[[7tQI(1 + f'{Ro))f + 3^}' 



(33) 



3 (3Z - 1) {[7rQ/(l + f(Ro))f + S^} - AirJ ^ VZ{Z - l){l + f'{Ro))Sl 

9l^Z(Z - 1) { [nQlil + r{Ro))f + S^} 

U = (^] = ^7T^QlJ[l + f{Ro)? 

\9QJsj 3^Z{Z - 1) {[rrQlil + /'(i?o))]' + S^} ' 



(34) 
(35) 



Numerical calculations show that the intensive quantities calculated by Eqs. ([33|) -([35 |) coincide with Eqs. ([24|) . (|25 
and pop , respectively. Thus, these thermodynamics quantities satisfy the first law of thermodynamics 



dM = TdS + VtdJ + UdQ. 



(36) 



IV. THERMAL STABILITY OF BLACK STRING 



Finally, we investigate the thermal stability of rotating black string solutions in f{R) gravity coupled to a matter 
field. The stability of a thermodynamic system with respect to small variations of the thermodynamic coordinates is 
usually performed by analyzing the behavior of the entropy S{M, J, Q) around the equilibrium. The local stability in 
any ensemble requires that S{M, J, Q) be a convex function of the extensive variables or its Legendre transformation 
must be a concave function of the intensive variables. The stability can also be studied by the behavior of the energy 
M{S, J, Q) which should be a convex function of its extensive variables. Thus, the local stability can in principle be 
carried out by finding the determinant of the Hessian matrix of M{S, J, Q) with respect to its extensive variables Xi, 
H^.j^. — [d'^M/dXidXj] [22,[23- I^i o^^ '^^-se the mass M is a function of entropy, angular momenta, and charge. The 
number of thermodynamic variables depends on the ensemble that is used. In the canonical ensemble, the charge and 
the angular momenta are fixed parameters, and therefore the positivity of the {d^M/dS^)j^Q is sufficient to ensure 
local stability. We find that the black string solutions are always thermally stable independent of the value of the 
parameters q and S. We have shown the behavior of {d^M/dS'^)j,Q as a function q and r+ for different value of 5 
and /'(i?o) in figures 4-6. These figures show that the black string solutions in /(i?)-Maxwcll theory with constant 
curvature scalar are always thermally stable. 
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FIG. 4: The function {d'^M /dS'^)j,Q versus g for i = 1, H = 1.25, r+ = 0.7 and Ro = -12. f'{Ro) = (bold line), f'{Ro) = 1 
(continuous line) and /'(i?o) ~ 2 (dashed line). 




FIG. 5: The function {d'^M/dS^)j,Q versus r+ for / = 1, H = 1.25, Ro = -12 and f'{Ro) = 1. g = 0.5 (bold line), q = 1 
(continuous line) and q = 1.5 (dashed line). 




FIG. 6: The function {d^M/dS^)j,Q versus g for i = 1, /'(i?o) 
line) and H — 2.25, (dashed line). 
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V. SOLUTION WITH NON CONSTANT RICCI SCALAR 



In this section we would like to extend the study to the case where the Ricci scalar is not a constant, instead we 
reconstruct it as i? = Ri^) as a result of our calculations. In this case, we find out that we can obtain solution only in 
the absence of the matter field. As we mentioned in the introduction, in general, the field equations of f{R) theory 
coupled to the matter field are very complicated and hence it is not easy to find exact analytical solutions. Thus we 
only consider the uncharged black string solution. We start with the following action 



1 
16^ 



d^x^^f{R) , 



where f{R) is a generic function of the Ricci scalar R. We also modify our metric (jH]) a bit as follow 



ds' 



-7V(r)e2"('-) (Sdi - ad(l>f + r^ (^ 



dt — So 






-dz^ 



(37) 



(38) 



where we have added an additional function a(r) in the metric coefficients. For simplicity we only consider the 
non-rotating black string with a ~ 0, thus the above metric reduces to 



ds' 



-7V(r)e2"('')dt2 



The scalar curvature for metric pQ]) reads 

dN [r) da (r) 



i? = -3- 



dr dr 



-2N{r) 



dr^ 
N{r) 



— a{r) 
dr 



-dz\ 



(39) 



d^N (r) 
dr"^ 



-2N{r) 



d^a [r) 
dr^ 



4 dN (r) AN (r) da (r) 2N (r) 
r dr r dr r^ 



(40) 



We use the Lagrangian multipliers method. In the fram.ework of Friedmann-Robertson- Walker universe this method 
was studied in |25l-l27|. while for static spherically symmetric black hole solutions it was investigated in [ij,[23- In 
this approach one may consider the scalar curvature R as independent Lagrangian coordinates in addition to the 
functions a{r) and N{r), which appear from the metric line element. 

Introducing the Lagrangian multipliers A, after using PO]) . the action p7p can be written 



S 



1 
16^ 



+2N{r) 



dt I dr 
d 



(e"('^V2)|/(i?)-A|i? + 



dr 



N{r) 



dr 



a (r) 



dr 



a (r) 



d^N (r) 



AN (r) da (r) 2N (r) 
r dr r^ 

Varying the above action with respect to R, one gets 



dr^ 
}}} 



2 7V(r) 



d^a (r) 4 dN (r) 



dr"^ 



dr 



A - fiR), 



(41) 



(42) 



where the prime denotes the derivative with respect to the scalar curvature R. Substituting this value and integrating 
by part, the Lagrangian takes the form 



L{a,da/dr,N,dN/dr,R,dR/dr) = 



[f{R) - RfiR)] - 2/'(i?) 



dN{r) 



dr 



N{r) 



dr \ dr 



2N{r) 



da{r) 
dr 



Making the variation with respect to a, one gets the first equation of motion 



Rf{R)-.f{R) , 2 
fiR) 

2N{r)f"{R) 



d^R 
dr^ 



N{r) 



dN{r) 
dr 



(43) 



(44) 



dN{r)/dr\ dR , f'"{R) 
d^ 



2N{r) 



fiR) 

The variation with respect to N{r) leads the second equation of motion 



fiR) 



dR 
dr 



0. 



'da{r) n"{R)dR\ 
dr \f'{R) dr) 


fiR) d^R 
f'{R) dr^ 


.r{R) fdR 
fiR) \dr 



= 0, 



(45) 



while by making the variation with respect to R, we recover Eq. (HUl) . Given f{R), together with equation (pn|. the 
above equations form a system of three differential equations for the three unknown quantities a{r),N{r) and R{r). 
We would like to note that one advantage of this approach is that a does not appear in Eq. pi| . In what follow, we 
will find exact solutions of the above system of differential equations. 

In the special case of constant curvature R = Rq and a = constant, it is easy to show that the only solution of Eqs. 
(|40|) and (j44|) is the Schwarzshild de Sitter black string solution with flat horizon. 



AT. N 2m A 2 

r 6 



(46) 



where m is a constant of integration which can be interpreted as the mass parameter of the black string and we have 
defined 2A = Rq — f{Ro)/f'{Ra) and Rq = 4A. Notice that in the absence of the matter field, q — 0, solution (fT2|) 
coincides with the result obtained in pS)) . as expected. 

Next, we consider the case of non constant Ricci curvature, but still with a = constant. From Eq. (j45p we have 



,,„ I dR 
•' ' dr 



-n''^) = ^m) 



(47) 



10 

which has the following solution, 

f{R)^mr + n, (48) 

where m and n are two integration constants. Given the explicit form of R, we may find r as a function of Ricci scalar 
and reconstruct f'{R) realizing such solution. From Eq. (j40l) with constant a, one gets 

^ ^ d^Njr) _ J dNjr) _ 2N {r) ^ ^^^^ 

dr"^ r dr r^ 

Using the fact that {f"{R))dR/dr = df{R)/dr = m and df{R)/dr = f'{R)dR/dr, and multiplying Eq. g^ by f'{R), 
we arrive at 

d^N{r) r n\ 4mN{r) 2nN{r) mdN{r)_ 

dr^ \ r J r^ r^ r dr 

When m — 0, the solution of the above equation is ones obtained in (j46p . For n — 0, the general solution is 

N{r) = -Cir2 + ^. (51) 

Substituting in Eq. (HHl) we again arrive at constant Ricci scalar, R — 12Ci. Although in this case f'{R) = df{R)/dR = 
mr is not a constant, but still we have df{R)/dr — 0, which implies that f{R) = constant. Next we look for the most 
general solution of Eq. ([5n)) with n 7^ and m 7^ 0. Solving ([501 . we find 

iV (r) = -Cir^ + — [2 n^ - imn^r + 6 mV^n - GmV^ In fm +-) 1 , (52) 

where Ci and C2 are two arbitrary constants. Given solution (|52|) one can basically construct f{R) by using Eqs. 
(|^5)) and (P^j) . In order to simplify the above solution, we choose n — \ and C2 = — l/m, 

Af(r) =3 6mr-Cir2 + 6mV2ln[m+- ) . (53) 

mr \ r J 

For this general case the Ricci scalar becomes 

i?(r) = 12Ci-72m2 In m+-+^ — >-. (54) 

\ r J r'^(mr + ly 

which is clearly not a constant. Now we want to reconstruct the corresponding f{R) theory. From Eq.(j48)), for n = 1 
one has 

Integrating (|55l) . by using ([54l) . we get 

„ / 1 \ 3 (I2m3r3 + ISm^r^ + 4mr - 2) 
/[i?(r)] --36m2 In m+-+^ 2^—7^2 ^ (^6) 

Combining Eqs. ([5^ . ([SS]) and ([5B|) . one gets the following differential equation for function f{R), 

/(i?)-#+6Ci = 0. (57) 



[/'(i?)]2[/'(i?)-l]2 ■ ^^ ^ 2 

This equation has a simple solution as 

f{R) ^ I - 6C1 - 48m2, (58) 

but it has also another complicated solution which we have not presented it here. Thus we have found a black string 
solution in f{R) gravity with non constant Ricci scalar. This approach also leads to construct Ricci scalar as a 
function of r, as given in Eq. ([M)) . The obtained solutions in this section differ from that presented in ^8] for axially 
symmetric solutions in f{R) gravity. It is worth mentioning that metric p8l) has a good property for which its static 
and rotating solutions coincide and so in this section we only study the static case. Following the approach of this 
section, one can easily check that solution (|52]) can be deduced for rotating case where a 7^ 0. Besides, in this section 
we only considered the case with a = constant, and derived the metric function (|52p as well as R{r). The study can 
also be generalized to the case where a = a{r). We leave it and also thermodynamic considerations of the obtained 
solution in this section, for future investigations. 
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VI. CONCLUSIONS 

In order to obtain the constant curvature black hole solution in f{R) gravity theory coupled to a matter field, the 
trace of the energy-momentum tensor of the matter field should be zero |14| . Since the energy- momentum tensor of 
Maxwell field is traceless in four dimensions, therefore spherically symmetric black hole solutions from f{R) theory 
coupled to Maxwell field was derived in four dimensional spacetime 14] . 

In this paper we continued the study by constructing a new class of charged rotating solutions in /(i?)-Maxwell 
theory with constant curvature scalar. This class of solutions describe the four dimensional charged rotation black 
string with cylindrical or toroidal horizons with zero curvature boundary. These solutions are similar to asymptotically 
AdS black string of Einstein-Maxwell gravity with suitably replacement of the parameters. However, the solution 
presented in this paper has at least two differences from AdS black string solutions of Einstein-Maxwell gravity. First, 
the conserved and thermodynamic quantities computed here depend on function f'{Ro) and differ completely from 
those of Einstein theory in AdS spaces. Clearly the presence of the general function /'(i?o) changes the physical values 
of conserved and thermodynamic quantities. Second, unlike Einstein gravity, the entropy does not obey the area law 
for black string solutions in /(i?)-Maxwell theory as one can see from Eq (f^ . We studied the physical properties 
of the solutions and found a suitable conterterm which removes the divergence of the action. We obtained mass and 
angular momenta of the string through the use of conterterm method. We also derived the entropy of the black string 
in f{R) gravity which has a modification from the area law. We obtained a Smarr-type formula for the mass namely 
M{S, J, Q) and checked that the obtained conserved and thermodynamic quantities satisfy the first law of black hole 
thermodynamics. Finally, we explored the thermal stability of the solutions in the canonical ensemble and showed 
that the black strings derived from f{R)- Maxwell theory are always thermally stable. This is commensurate with 
the fact that there is no Hawking-Page phase transition for black objects with zero curvature horizon [24] . 

We also extend the study to the case where the Ricci scalar is not a constant. For this purpose we used the 
Lagrangian multipliers method and found an exact black string solution in f{R) gravity. In this approach one may 
consider the scalar curvature R{r) as an independent Lagrangian coordinates in addition to the metric functions and 
deduce R{r) as a solution of the field equations. We found the explicit form of R{r) as well as the metric function 
which has a logaritmic term. It is worth noting that since for the non constant Ricci scalar, the field equations of f{R) 
theory coupled to the matter field become very complicated, in this case, we could only derived analytical solution in 
the absence of the matter field. 
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